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O ■ Abstract 

■ 

• • . We calculate the rate of energy emission by bremsstrahlung of neutrino pairs 

^ . by electrons moving in the crystalline lattice of ions in dense matter in the 

crust of a neutron star. Since the periodic potential in the solid gives rise to 
electronic band gaps which can be as large as about 1 MeV, it is necessary 
in estimating the bremsstrahlung rate at low temperatures to take into ac- 
count band structure in detail. We find that, in the densest parts of the inner 
crust of a neutron star, neutrino emission at temperatures of about 2 x lO^K 
or less is much suppressed compared with earlier estimates that treated the 
electron-lattice interaction perturbatively, and conclude that neutrino pair 
bremsstrahlung by electrons in the crusts of neutron stars is much less impor- 
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tant for neutron star thermal evolution than was previously thought. 
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I. INTRODUCTION 



The possible importance of neutrino pairs generated in collisions of electrons with nuclei 
was pointed out by Pontecorvo , and the first detailed calculations of the rate of the pro- 
cess were made by Gandel'man and Pineev [Q], who considered the case of non-degenerate, 
non-relativistic electrons scattering from uncorrelated nuclei. This work was extended by 
Festa and Ruderman 0, who performed calculations for degenerate electrons. They al- 
lowed for relativistic effects, and took into account correlations between ions and screening 
of the ionic Coulomb potential by electrons. Their results were extended by Cazzola, De 
Zotti and Saggion 0], who performed calculations for uncorrelated ions for a range of con- 
ditions. A general many-body formulation of the problem of emission of neutrino pairs by 
electrons scattering from ions was made by Flowers ||^, who applied his results to estimate 
bremsstrahlung by scattering from the static lattice and also gave expressions for rates of 
processes in which lattice phonons are emitted or absorbed. Flowers also stressed the im- 
portance of the finite nuclear size in reducing the bremsstrahlung rate. The effects of weak 
neutral currents was investigated by Dicus, Kolb, Schramm, and Tubbs 0. Further stud- 
ies were made by Soyeur and Brown 0, who compared energy loss by the bremsstrahlung 
process with other energy loss mechanisms for various neutron star models. In recent years 
the process has been studied in detail by Itoh and collaborators |^. 

Basic to all the above works is the assumption that the effects of the electron-ion inter- 
action may be treated in lowest-order perturbation theory, and the two perturbation theory 
diagrams considered are shown in Fig. 1. Since the energy and momentum carried away by 
the neutrino pair are of order ksT and ksT /c respectively, and since the final particles are 
on mass shell, it is easy to see that the intermediate energy denominator must be of order 
ksT. The dimensionless parameter in the perturbation expansion is thus of order Vq/ksT, 
where Vq is the Fourier transform of the potential due to the ions. For a random collection 

of ions, the Fourier transform of the potential experienced by an electron will be of order 
1 1 

n^Vq/fl^ where Vq is the Fourier transform of the potential due to a single ion, nz is the 
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number density of ions, and fl is the volume of the system. Thus one does not expect prob- 
lems with the perturbation expansion in this case. However, the situation is very different 
if matter is a solid, with ions ordered on a periodic lattice. The Fourier transform of the 
electron- ion interaction is then of order uzVk, where K is a reciprocal lattice vector and 
this will lead to energy splittings in the electron spectrum of this order of magnitude. When 
the temperature is low compared with this splitting, perturbation theory fails. As we shall 
argue in detail below, the splittings can be as large as 1 MeV in the inner crust of neutron 
stars, and therefore at temperatures below about 10^° K the first-order perturbation theory 
calculations of the rate fail. 

In this paper we calculate bremsstrahlung of neutrino pairs by electrons moving in the 
periodic potential of the ions. Our primary purpose in this paper is to examine elastic 
scattering of electrons from the lattice, since this process was previously thought to be 
more important at low temperatures than processes in which lattice phonons are emitted or 
absorbed Instead of using perturbation theory to treat the electron- ion interaction, we 
first calculate band electron states, thereby implicitly taking into account the electron-ion 
interaction to all orders in perturbation theory, and then evaluate the rate of emission of 
neutrino pairs when an electron makes a transition from one band state to another. What we 
find is that bremsstrahlung of neutrino pairs is suppressed exponentially at temperatures low 
compared with the strength of the ionic potential, which determines the splitting between 
bands. A brief account of part of this work has been given previously 

The paper is organized as follows. In Sec. II we describe the basic physics and derive 
an expression for the neutrino bremsstrahlung emissivity. Analytical expressions for the 
emissivity in the high- and low-temperature limits are derived and numerical calculations 
are made for intermediate temperatures in Sec. III. Section IV contains a discussion of the 
results, and a comparison of neutrino emission by bremsstrahlung from the periodic lattice 
with that from conversion of phonons into neutrino pairs. 
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II. BASIC CALCULATIONS 



The basic process we consider is shown in Fig. 2, where the double hnes represent a 
band electron. If the band electron propagators in Fig. 2 are expanded to first order in the 
electron- ion interaction, one recovers the two diagrams shown in Fig. 1. To lowest order in 
the weak interaction, the rate for an electron in state 1 to emit a neutrino pair in a transition 
to a state 2, assumed to be unoccupied initially, may be found from Fermi's Golden Rule, 
and is given by 

dn='^6{Ef-E,)\Hf,\^dT, (1) 

where Hji is the matrix element of the weak interaction between the initial electron state i 
and the final state /. The energies of all particles in the initial and final states are denoted 
by Ei and Ef. We denote the energies of the initial and final electrons by Ei and E2 
respectively, and those of the outgoing neutrino and antineutrino by E^, and Ep, and dT is 
an element of the volume of phase space. The total rate of energy loss in neutrino pairs per 
unit volume is thus given by 

E = §^j:m- E{)m){i - mm^^wE, + e,). (2) 

Here f{E) is the Fermi function, and the factors f{Ei) and 1 — /(-E2) in Eq- (@) take into 
account the Pauli principle, since the transition can occur only if the initial electron state is 
occupied, and the final one empty. There are no blocking factors for the neutrinos, since we 
assume that neutrinos are able to escape easily, a good assumption for neutron stars except 
for some seconds immediately following the birth of the star. 

To calculate the rate at which an electron in one band state makes a transition to another 
band state we need to evaluate matrix elements of the weak interaction Lagrangian 

£ = -^^G^.Tl^.^iCLTl + Cj,T%)^, (3) 

using band states for the electrons, rather than the plane waves of most interest in the 
majority of apphcations. In Eq.(^), G is the Fermi coupling constant, = 7"(1 =F75)/2, 
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and in terms of the weak mixing angle 9w, Cl = 1 + 2 sin^ 9w and Cr = 2 sin^ 9w for the 
emission of electron neutrinos. The corresponding couplings for the emission of muon and 
r neutrinos, are = — 1 + 2 sin^ 9w and = 2 sin^ 9w- 

For the moment, let us be quite general, and consider transitions in which the initial 
electron state is specified by the wavenumber, p, a band index, n, and a spin index, a, 
and the final state is specified by similar quantities with primes. The squared modulus 
of the matrix element, summed over spins of the initial and final electron, neutrino and 
antineutrino states, is given by 

H = i2G27^a/3(pl,(Tl,nl|Jq^+qJp2,a2,n2)(p2,a2,n2|J^q,_qJpl,(Tl,nl) , (4) 



where 

7" — _ 

and the neutrino part of the trace is given by 



J^ = ^J d'xe-^'^-iCi^ri + CnT^) , (5) 



= E u^A^i)riv^Mv,M^lu^M . (6) 



The neutrino part may be evaluated by using the fact that for ultrarelativistic fermions 

u 

u,,{k)®u^^{k) = ^, (7) 

and one finds 

= 771^(^92^ + Q1Q2 - g'^^qi ■ 12 + *e"^^^(gi),(g2)5) • (8) 
lqiiq2 

The electronic matrix elements depend on the neutrino momenta only through the total four 
momentum of the pair, q = (tu/c, q) = qi + q2 and therefore one may sum over neutrino states 
in the expression for the energy emission rate, Eq.(0), keeping the total four momentum of 



the pair fixed. We then use the result |10 



E -E.- = ^(2g"/ + 9'^^')9iu - c|q|) , (9) 



to show that 

E ^q,qi+q.^(^ " " ^p)^"'' = " " c|q|) . (lO) 

qi,q2 ^'^^ 

Here 6'(x) is the Heaviside step function, which is unity for x > and vanishes otherwise. The 
form of this result, apart from the numerical coefficient, follows directly from the observation 
that the left hand side of this equation is a second rank Lorentz tensor with vanishing 
divergence. The vanishing of the divergence is a consequence of neutrino current conservation 
for free neutrinos, as may be seen by contracting Eq.(6) with q". Inserting Eq. (p!oD into 
Eq.(l), we find 



Pij P2 

((l|g-Jq|2)(2|g-J_q|l)-g2(l|J^"|2)(2|J_qJl)). (11) 



This expression is quite general: effects of band structure are included, through the matrix 
elements, and no assumption has been made regarding how relativistic the electrons are. 

First we make a general remark about what transitions are possible. The bremsstrahlung 
process is kinematically forbidden for initial and final electron states in the same band, n. 
The argument is essentially identical to the one to demonstrate that (photon) bremsstrahlung 
by a free electron is impossible. As a consequence of the fact that the electron (group) 
velocity VpE cannot exceed c, the energy difference, Ei — E2 between electron states must 
be less than cq, where q is the total momentum of the neutrino pair. However, this energy 
difference, which is equal to the energy of the neutrino pair, must exceed cq. Consequently, 
it is impossible simultaneously to conserve energy and momentum for electron states in 
the same band. This argument does not apply to transitions between electron states in 
different bands because the energy difference between the two electron states will generally 
be non-zero for small momentum transfers. 

We now turn to the case of the inner crust of neutron stars, where the electrons are 
highly relativistic, with Lorentz factors of 50 or more. Under these circumstances it is an 
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excellent approximation to neglect the mass of the electron, since this leads to errors of order 
(meC^//ie)^, which is less than 10^'^. Here fie is the electron chemical potential (neglecting 
electrostatic contributions), which in the extremely relativistic limit is simply ppc, with pp 
being the electron Fermi momentum. The electron helicity is a good quantum number in 
this limit, and this leads to several simplifications. First, the left- and right-handed currents 
satisfy a continuity equation g ■ Jq = for each helicity, and therefore the contributions in 
Eq.(|lT]) involving such quantities vanish. The energy emission rate is thus 

G2 



P15 P2 

6{E, -E2-0J) f{E,) (1 - /(E2)) 9{u; - c|q|). (12) 



Second, in the term of the form J ■ J the two hehcities decouple, and the net result is 
proportional to Cl + = 1(Cl + C^), where Cy = {Cl + Cr)/2 and Ca = {Cl - Cr)/2. 
The rate of energy emission is thus (| + sin^ Ow)'^ times the rate calculated in V-A theory, 
for which Ca = Cy = 1. The rate in muon and r neutrinos is obtained in a similar fashion. 
For them the vector and axial vector couplings are C'a = 1 — Ca and Cy = 1 — Cy [^] . The 
total rate of neutrino energy emission in all types of neutrinos is thus given by replacing 
C2 + Cl with Cl + C2 + 2((1 - Ca^ + (1 - Cyf). 

The next task is to evaluate the weak interaction matrix elements, and to do this we 
require a knowledge of the electronic states. Since splittings between bands are small com- 
pared with a typical electron energy, it is a very good approximation for all but a few electron 
states to use the almost-free-electron approximation, in which one takes into account just 
one Fourier component of the periodic potential at a time. The component of the lattice 
potential with wavenumber K mixes free electron states with wavenumber k with those 
with wavenumbers k±K. In the nearly-free electron approximation (see e.g., Ashcroft and 
Mermin [rT|), one considers only the mixing of the two electron states that are closest in 



energy, and for definiteness we denote the wavenumbers of the two states by k and k — K. 
The energy eigenvalues are given by 
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E^ik) = ± y'(!!LZiiz5i). + v^ , (13) 

where the plus sign corresponds to the upper band, and the minus sign to the lower band. 
Here ek is the energy of a free electron, and Vk is the matrix element of the static part of 
the electron-ion interaction between the two free-electron states that are mixed. Since the 
helicity of the electrons is conserved in the ultrarelativistic limit that we arc treating, the 
matrix element is simply the Fourier transform of the potential for wavenumber K acting on 
an electron, multiplied by a factor giving the overlap between the states with wavenumbers 
k and k — K, but of the same helicity. The helicity factor is, apart from a phase, just 
cos(|^k,k-K), where ^k,k-K is the angle between the two electron wave vectors. In our 
calculations the band states of greatest interest are those close to the Fermi surface, and 
with a component of the wavevector along K close to ±K/2. For this case the modulus of 
the helicity factor is then simply v± — {1 — {K / {2kF)Y)^ . The Fourier component of the 
electrostatic potential due to the ions is given by Ane / {p\^) , where (pk); is the statistical 
average of the Fourier transform of the charge density of the ions. As a consequence of 
screening by electrons, the total electrostatic potential due to both ions and electrons is given 
by the expression for ions alone, reduced by a factor of the electronic dielectric function, 
6k- In the standard Thomas- Fermi approximation this is given by £k = 1 + {kxp/KY, 
where for the screening wavenumber we use the result for an ultrarelativistic electron gas, 
k'^p — Aakp/ir, with a — /{he). In summary, the matrix elements of the lattice potential 
are given by 

where riz and Z are, respectively, the number density and proton number of the ions, Fk 
is the form factor reflecting the shape of the nuclear charge distribution and and W is the 
Debye- Waller factor describing thermal vibrations. 

The evaluation of the electron trace is more complicated because of the electron band 
structure. First of all, a remark about bookkeeping: the electron spectrum is shown in Fig. 4 
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as a function of k\\ the component of the electron wavenumber parallel to K for a fixed value 
of k±, the component of the electron wavenumber perpendicular to K. For simplicity we 
take into account only the effects of the lattice potential at a single reciprocal lattice vector. 
At low temperatures the electron states that can participate in neutrino pair emission are 
those close to the maxima and minima of the electron spectrum. With the normal choice of 
Brillouin zone boundaries, these states would lie close to the zone boundary, and one could, 
for example, have processes in which an electron on branch 1 made a transition to a state 
on branch 2, or to a state on branch 4. The fact that the wavenumbers of states on branches 
2 and 4 differ by K complicates the calculations, and it is convenient to redefine the zone 
boundary so that the branches 3' and 4' occur in the same zone as branches 1 and 2. All 
states with k\i ~ K/2 are then given by a single expression, with no discontinuities. For 
definiteness the wavenumbers of the two plane wave states which are strongly mixed will be 
denoted by k and k — K. 

The expression for the electronic part of the matrix element may be simplified by making 
use of the fact that we are interested in processes at low temperatures, and therefore the 
energy and momentum of the neutrino pair are small compared with the Fermi momentum 
of the electrons, and the reciprocal lattice vectors. In particular, the small momentum of 
the neutrino pair means that the only Fourier transforms of the current operators of interest 
are ones with long wavelengths, and therefore q = Pi — P2- Matrix elements of the current 
operator between the two states will possess components with wavenumbers pi — p2 ± hK., 
Umklapp terms, but these will not enter in our calculation as a consequence of the fact that 
the low temperature restricts possible neutrino momenta to values much less than hK. We 
note that had we made the conventional choice of position of the zone boundary, we would 
have been forced to consider processes in which pi — p2 is comparable in magnitude to hK. 
This reflects the fact that what are normal processes for one choice of the zone boundary 
may be Umklapp processes for some other choice. The states in the upper and lower electron 
bands are 



10 



M/p,,(r) = -l=(v^e^nP'^uM - «pe^(P-^^>^«.(p - ^K)) . (15) 



The expressions for the "coherence factors" Up and Vp will be given below. 

In calculating long- wavelength components of the current matrix elements it is convenient 
to invoke the identity 

Mai(k)7^M^2(k) = — 5<,,,^2 ■ (16) 

We find 

(P2,cr2,n2|jq"|pi,a-i,ni) 



^q„pi — P2 ^ni ,n2+l — ,0"1,(T2 + Cr5+ 

,ai,cT2) ~'r Umklapp , (17) 



Pf 

where we have neglected terms of order Ap/pp, where Ap is the difference between the 
electron (pseudo)momentum and the Fermi momentum. In Eq.(0), Mi 2 = Mpi a; "^1,2 = 'i^pi,2 
and the electron four momenta pi^2 = (-£'i,2/c, Pi,2) occur in the combinations P = {pi + 
^2)72 , P' = P — {0, hK.). We thus find for the electron trace 

(1| Jq|2)(2| J_q^|l) = ^-^ ^Mlt^ltt2^^2^q,pi-p2^ni,n2 + l + Umklapp , (18) 



where v\\ = yl — f ^ = K/(2kp). From Eq.(^) follows that the energy emission rate in 
neutrinos is given by 

G2 ry2 i „ 

^ 7f,io ^ 2^ ^11 / 9d Pid P2 ujq U1V2V1U2 

5{E, -E2- uj)5^^\vi - P2 - q)/(^i)(l - f{E2))e{uj - c|q|) . (19) 

The sum over electron bands is contained in the sum over reciprocal lattice vectors K in 
Eg. (p!9D . Note that since electron states were defined within a given Brillouin zone, Eq.(|19|) 
includes a factor of 1/2 so that there is no overcounting when summing over equal but 
opposite reciprocal lattice vectors. From Eq.(p!7|) one also finds that (l|g- Jq|2) (2|g- J_q|l) cx 
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{uiV2q ■ P — viU2q ■ P'Y . It may be verified that this vanishes identically, as argued earlier 
on the basis of current conservation. 

We now perform the energy and momentum integrations. To do this we require expres- 
sions for the "coherence factors", Up and fp. They are identical in form to those found in the 
standard treatment of nonrelativistic electrons in a weak periodic potential and are given 
by 



K 

Up = I = , ^^p 



< = kl + vl = 1(1 - ^), and up^p = ^ , (20) 



with ^p = (ep — ep_;iK)/2 and £p = + V^. We also note that, in terms of the variables 
(P^), the energies of the upper and lower electron states are 



E^{Ap^, Ap") = pfc + Ap^v^c ± ^J{Aph\\cy + , (21) 

from which we find u = cvj_q± + Si +82- 

It is convenient to express momenta in cylindrical coordinates, with respect to K. Also, 
since the electrons participating in neutrino pair emission are close the maxima and minima 
of the electron spectrum, we define the shifted variables 

Ap" = p" — v\\Pf and Ap"*" = p"*" — v±pf ■ (22) 

We also have 

q±=pi -P2 = ^Pi - ^P2 

Q\\ =Pi-P2 = ^Pi - ^pI 

q<l> = PFVj_{(t>l - 02) , (23) 

where in the last expression, we have discarded curvature terms, which contribute only to 
next-to-leading order in ksT /ppc. The p2 integration is now trivially performed, with the 
aid of momentum conserving delta function, and the pi angular integration simply gives the 
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constant 27r. Since the coherence factors do not depend on p^, the integration over is 
straightforward and yields the result 

J{E,-u;)-f{E^) 



J pidpif{E,){l - f{E,)) = v^pF j dApi 



-"--c—l • (24) 

where we have again utilized the fact that characteristic temperatures are much smaller than 
the Fermi energy. Performing the remaining energy integration finally yields the expression 
for the energy emission rate per unit volume, 

^ ^ 247r6^^°c3 ^2 ^ ^^^^ U1V1U2V2 _ ^ 9{u - c|q|) , (25) 

with u! = El — £'2- 

Inserting into Eq.(P3|) the expressions for the coherence factors (^Dj) and using Eq.(PTD 
we find the energy emission rate, 

with P|| = Ap\ — ^q\\is an average crystal momentum variable for the initial and final 
electron states. The integral in Eq. (|26|) may now be performed, leaving an integral over 
the remaining momenta, which we express in terms of the dimensionless variables 

c^ll cgii cq± 

The energy emission rate may thus be written as 

E = J^I^^^A±^^^^kBTfy^vlV^I{v.,'^) , (28) 
567/i^°c9 2 « ; II K V ll'|y^|^ ' ^ ^ 

where J is the dimensionless integral 

^M = ^ d^i / / . c^X3^ . J ■ 29 

Stt' jo Jo Jx™" 6162 (e'" — 1) 

The integral is to be carried out over the region where — — is positive. We have 
defined the dimensionless energies 
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ei,2 = Y ^ + vf\{xi ± ^^2)^ and w = v±X3 + 61 + 62 , (30) 

in terms of the variable t = ksT /V\^. The hmits of the X3 integration in Eq.(P5D are 

max,min ^-L / , \ i IT . \9 2 2 /oi\ 

= -2(61 +62) ± ;^ y/(6i +62)^ . (31) 

The hmits may be put in a more symmetric form by defining the variable y = — v±{ei + 
62)/^^, leading to 

/= — -vfj dxi dx2 / dy(a^-y^)2 ^ , (32) 

StT^ II Jo Jo 6162 J-a ^ ' v^y+^{e^+e2) ' ^ ^ 

e II — 1 



with a = jijiei + 62)^ - vf^xl 



III. EVALUATION OF THE NEUTRINO EMISSION RATE 

In this section we evaluate our basic result for a number of different situations, beginning 
with the low- and high-temperature limits, and then we describe results for more general 
cases. 



A. High-and Low- Temperature Limits 

Analytical results may be obtained in limiting cases. We first consider temperatures 
small compared with |Vk|, where our results differ dramatically from earlier ones. This 
limit corresponds to t <^ 1. In this case, the denominator in Eq.(^U|)is essentially e"^, and 
we simply expand the integrand around the state of lowest "energy" w, given by 

xi = X2 = and X3 = xf = - ^ ^ j ■ (33) 

All three integrals are thus reduced to gaussian ones, with the result 

/ = ^t-2e^3 . (34) 

The rate of energy emission is therefore 
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The exponential dependence reflects the fact that the minimum energy of the neutrino pair 
is 2|Vk|/(1 + v±). This is easily seen by observing that the energy of the neutrino and 
antineutrino is given by = Ei — E2 > 2|Vk| + <l±v±, and that, in addition, the four- 
momentum of the neutrino pair must be time- like, uj > cq > c\q±\. 

We next consider the limit of temperatures much greater than Vk, corresponding to 
t ^ 1. In this limit, the range of y integration in Eg . (1321) goes to zero in the region 
2xi < X2- It is furthermore convenient to perform the Xi integration last, since for fixed xi, 
X2 and y lie on a disk of radius 2xi. We write 



27r7 Jo J^2x, - 1 

where 

1,= !^' i,, i^A-^-v 



^ > 2 



^0 Axf — X2 

= ^(4x?-3y^ + M!) . (37) 

The remaining integrals are performed by defining the variables rji = 2xi + v±y, associated 
with energy, and the orthogonal variable 772 = —'v±xi + 2y. The 7]2 integral is elementary, 
and the final, rji, integral is evaluated using standard techniques for Bose integrals 



00 in^ 







dm — '-^ — = . (38) 

' e'^i - 1 63 ^ ^ 



We are thus led to the following result for the basic phase space integral 



1 H , 1 



1-:! log 7:2) ' (39) 



2 V 2 o 2 

V,V\\ V, V 



with the corresponding energy emission rate given by 



We have checked that Eq.(^) is consistent with earlier calculations in which the electron- 
lattice interaction was treated perturbatively ( Fig.l ). Flowers has given the energy 
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emission in integral form, valid for nonvanishing electron mass. By evaluating these integrals 
in the limit of vanishing electron mass we were able to arrive at the analytical result ( |^ ) . 



B. Results for Arbitrary Temperatures 

The general expression for the neutrino emission rate is rather complicated, so we begin 
by discussing a number of its ingredients. 

In our discussions we shall focus on the highest density regions in a neutron star crust 
since most of the mass of the crust is at densities close to that at the inner boundary of the 
crust, and consequently such matter is expected to be responsible for most of the neutrino 
pairs produced by electron bremsstrahlung. An important quantity is the strength of the 
periodic potential, given by Eq.(|l^. First of all, we note that the effects of electron screening 
are negligible, since for the smallest reciprocal lattice vector for a bcc lattice, K = 2\^n/a, 
and therefore the dielectric function is e = 1 + {9/{27T^)y/^Z'^/^a ^ 1 + 1.79 x lO^^Z^/^, 
which differs from unity by less than 3% for the largest Z's encountered in neutron star 
crusts. Next, we turn to the nuclear form factors. At the highest densities, their effects are 
significant, as we now show. According to the calculations of Ref. [|l^, spherical nuclei have 



a radius of r^v = 8.7 fm when spherical nuclei become unstable with respect to formation of 
rod-like nuclei. The interface between nuclear matter and the neutron liquid outside nuclei 
is rather diffuse, so we assume that the nuclear charge distribution may be approximated 



by a gaussian exp(— 3r^/(2r^^g)), and for Vrms we take the value yS/Sriv, corresponding to 
a liquid drop with uniform proton density. The charge density in coordinate space is thus 
proportional to exp(— 5r^/(2r^)), and the form factor is given by Fg = ea;p(— g^r^/10). For 
matter at the highest density at which nuclei are roughly spherical, for which /le = 78MeV 
and Z = 48 [|12|[| , the form factor is 0.675 for the smallest reciprocal lattice vector, and 



^We draw attention to the fact that in Ref. [12| the cell radius at the phase transition between 
spherical nuclei and rod-shaped ones the cell radius is given as 19.2 fm. This is a misprint, and the 
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0.009 for the largest one, approximately 2pp. 

Finally, we consider the Debye- Waller factor. This is more difficult to estimate since, to 
the best of our knowledge, there exists no treatment of lattice vibrations in dense matter 
that takes into account the effects of the neutrons between nuclei. To make estimates we 
shall use the results for the pure Coulomb lattice given by Yakovlev and Kaminker W5 



Again we assume conditions appropriate for the highest density at which spherical nuclei 
survive. The total number of nucleons per unit cell is A' ^ 1466, and the number of nucleons 
within the nucleus is A ^ 427. Thus the total number of nucleons per proton, x is 0.033. 
(Note that, because there are neutrons outside nuclei, x ^ Z /A.) The total mass density is 
1.07 X lO^^g cm~^. The Debye- Waller factor is 

W{q) ^ 0.078-(— if fl.4e-9-" + 12.995t) {^f . (41) 

Here t = T/Tp, where the "plasma" temperature is given by 

^ kB ke^ rriion 0.1 0.IA/ ' ^ ^ 

where Up is the plasma frequency for the ions alone, neglecting the presence of the outside 
neutrons. The quantity pu denotes the mass density in units of 10^^ g cm~^. For a tem- 
perature of 10^ K, e~^^'^^ is 0.98 for the smallest reciprocal lattice vector, and 0.81 for the 
largest one. At lO^*' K, the corresponding values are 0.83 and 0.14 . This demonstrates that 
the effects of Debye- Waller factors can be significant for temperatures close to the melting 
temperature, which is about 1.2 x lO^'' K, but the effects of the periodic potential are still 
substantial for the lowest reciprocal lattice vectors. We emphasize that these estimates are 
based on the assumption that the lattice dynamics is that of a pure Coulomb crystal, and 
more realistic estimates can be made only after a thorough investigation of lattice dynamics 
that allows for the neutrons between nuclei. 

For conditions under which nuclear form factors and Debye- Waller factors may be ne- 
glected, the dependence of emission rate on the pe, T, and Z is of the general form 



actual value is 17.6 fm. In our earlier work S we used the spurious value, which led to a Z of 62. 
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i?~TV^^7(Z,T//ie) , (43) 

where g approaches a constant for temperatures large compared with Z^/^a/Xe, the magnitude 
of the largest band gaps. 

Before presenting the results of detailed calculations we describe qualitative features of 
the emission rate. As one can see from the asymptotic low- and high-temperature expressions 
Eqs.([35|) and (^OD , the crossover between the two sorts of behavior occurs at a temperature 
~ 2|Vk|/(1 + v±)/kB- Thus the crossover temperature for large wavevectors,~ 2pi?, is a 
factor ~ Z'^/^ smaller than for the smallest reciprocal lattice vectors, and thus there is a large 
temperature range in which neither of the limiting expressions is a good approximation. We 
note that reciprocal lattice vectors close to 2pF give very small contributions because helicity 
conservation ensures that electrons cannot be scattered by ions for K = 2pp, as reflected by 
the v± factor in the expression for the ionic potential, Eq. (|T^) . At high temperatures, a small 
reciprocal lattice vector gives a larger contribution than a large one, because the emission 
rate is proportional to the square of the Coulomb matrix element. As the temperature 
is lowered, the first reciprocal lattice vectors for which bremsstrahlung is suppressed are 
the smallest ones, because they give rise to the largest band gaps. Consequently, at lower 
temperatures, the reciprocal lattice vectors that give the largest contributions become larger 
and larger, until in the low-temperature limit, the dominant ones have wavevectors close to 

We now present results of a numerical evaluation of the the rate of energy emission from 
neutrino pair bremsstrahlung. We consider matter at the highest density at which nuclei are 
roughly spherical, for which /ig = 78MeV and Z = 48, according to the recent calculations 
of Ref . |12[ . We shall neglect Debye- Waller factors. 



In Fig. 5 we plot energy emission rates, for several values of v\\ = K/{2kp), compared 
with the corresponding quantities in the high temperature limit. The solid lines in Fig. 5 
are for nuclear form factor Fq = 1, and the dashed lines are for Fg = exp(— where 
r^r = 8.7fm. 
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The total energy emission rate is obtained by evaluating the basic phase space integral, 
Eq.(p2D for all reciprocal lattice vectors entering in the sum of Eq.(p6D. The relevant recipro- 
cal lattice vectors are those of length K < 2kp, where kp is the electron Fermi momentum. 
The Fermi sphere is of radius R = {3Z / {32tt)Y^^ in units of Arc /a, the side length of the 
conventional cubic cell of the reciprocal ( fee ) lattice. Here, a = is the length of a side 
in the conventional cubic cell of the direct ( bcc ) lattice. For Z = 48 we find that R = 1.13, 
and that there are 200 reciprocal lattice vectors contributing to the sum. 

The numerical evaluation may be streamlined by exploiting the symmetry of the lattice 
when evaluating phase space integrals. This becomes especially relevant when there is a 
large number of reciprocal lattice vectors contributing to the energy emission rate. A given 
reciprocal lattice vector has 47 partners of equal length, if the three vector components differ 
in magnitude, and none are zero. The degeneracy is reduced if two or more components are 
equal, or if one or more components are zero. To give an idea what is to be gained from 
such considerations, we note that the 200 reciprocal lattice vectors contributing to neutrino 
emission have only 11 distinct lengths. 

The three dimensional integral (^2]) was evaluated using the extended trapezoidal rule 

. We checked that both the low temperature ( [5^ ) and high temperature ( PD| ) limits were 
reproduced to within better than 1% for a range of values of v\\, corresponding to different 
choices for reciprocal lattice vectors. 

In Fig. 6 we show the rate of energy emission from neutrino pair bremsstrahlung, divided 
by the high temperature rate, Eq. (^OD , for the conditions we have described. The solid 
line shows this ratio for nuclear form factor Fg = 1. The short dashed line shows the 
corresponding ratio for Fg = exp(— g^r^/10), with = 8.7fm. The rate from the 12 
smallest reciprocal lattice vectors alone divided by the total high temperature rate is given 
by the line with alternating short and long dashes, for Fg = 1. The long dashed line gives 
the emission rate due to conversion of thermal phonons into neutrino pairs, Eq.(44), as 
compared with Eq.(40). 
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IV. DISCUSSION OF RESULTS 



An unexpected feature of the results is that the ratio E/T^ is not a monotonic function of 
temperature, and it has a maximum. We do not beheve this to be an artifact of the numerical 
calculations for two reasons. First, the numerical calculations agree with the analytical 
results. This is a non-trivial check since the regions of integration in the high- and low- 
temperature limits are very different. Second, we have performed the numerical integrations 
by two different methods, with identical results. It is interesting to note that non-monotonic 
temperature dependence was predicted theoretically and subsequently observed in a quite 
different physical context, the relaxation of nuclear spins in superconductors. We surmise 
that the origin of the non-monotonic behavior in the case of neutrino emission may be a 
similar effect, because the excitation spectrum of electrons near a band gap is similar to that 
for the excitations in a superconductor. However, the neutrino emission process we have 
considered is an interband process, and thus its analog in a superconductor would be the 
annihilation of two positive-energy excitations. However, the relaxation of a nuclear spin 
in a superconductor takes place by scattering of one positive energy excitation to another 
positive energy state, which is analogous to an intraband process in band theory. It is 
desirable to confirm this analytically, by, for example, calculating the leading corrections to 
the high-temperature limit. 

While the suppression of the bremsstrahlung rate for processes with a particular recipro- 
cal lattice vector is rather abrupt, the suppression of the total rate from all reciprocal lattice 
vectors is less abrupt, as a consequence of the fact that the characteristic temperatures for 
suppression span a range of order Z^/^. 

Prom Fig. 6 we see that the smallest reciprocal lattice vectors give almost 2/3 of the total 
rate at high temperature. As the temperature is lowered, this fraction drops rapidly, since 
the emission from the smallest reciprocal lattice vectors is the first to be suppressed. 

Figs. 5 and 6 also illustrate the effect of nontrivial nuclear form factors. For a given 
reciprocal lattice vector K, the crossover temperature is lower if the nuclear form factor 
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Fk differs from unity than if Fk = 1, since the matrix element of the lattice potential is 
smaller. For large where Fk can be very small, this effect is much more pronounced than 
for small K. However, the contribution of large reciprocal lattice vectors to the total rate 
is small, due to the v± factor in the ionic potential. Consequently, the magnitude of the 
change resulting from nontrivial nuclear form factors is closer to that given by the smallest 
reciprocal lattice vectors. 

In our calculations we have assumed that the periodic potential mixes only pairs of free 
electron states. However, for points of high symmetry in the Brillouin zone, it is possible for 
more that two free-electron states to be degenerate. At high temperatures such points will 
be of little significance because the phase space for transitions to states in their vicinities 
is small compared with the total. However, at low temperatures they could be significant 
if it happened that two bands could cross, even when the effects of the periodic potential 
are taken into account, since this would give the possibility of the neutrino emissivity at 
low temperatures behaving as a power of the temperature, rather than exponentially. With 
this in mind we have investigated the degeneracy of bands at points of high symmetry. The 
standard discussion of this problem for terrestrial solids cannot be carried over to case 



of dense matter because electrons are highly relativistic, and therefore spin-orbit coupling is 
strong. However, by generalizing the standard treatment to take into account electron spin 
and spin-orbit coupling by employing double groups, one can demonstrate that for the bcc 
structure, bands cannot be degenerate, and therefore the neutrino emissivity will behave 
exponentially at low temperatures. We conclude that points of high symmetry will not play 
an overwhelming role at low temperatures, and consequently our calculations that allowed 
for mixing of only two wavenumbers at a time will give a reliable estimate. 

Another process that can generate neutrino pairs in the crusts of neutron stars is the 
conversion of thermal phonons into pairs. This has been considered by many authors 
and most recently by Yakovlev and Kaminker |T^. One important characteristic tempera- 
ture is the "plasma" temperature defined in Eq.(H^), in the approximation that the effects of 



neutrons on the lattice dynamics are neglected, except insofar as they affect the relationship 
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between the total mass density and the electron density. Another important temperature is 
the melting temperature, = Z'^e'^/{r^kBT^) ^ 2.6 x lO^^iZ/GOf^^ipux/O.lY^^ K, where 
for the characteristic coupling parameter at melting, we have taken the value 172 ob- 
tained by Nagara et al. [jl^ The quantity Vc = (47rr7,^/3)~^/^ = (87r/3)"'^/'^a, where the latter 
result is for a bcc lattice. For temperatures larger than about Tp/16, the energy emission is 
given approximately by 

Ephonon ^ 4.2 X 10^(x/0.1)2/Vr4 ^'rjerg g-^ s''. (44) 

Thus at the melting temperature, the neutrino emission due to the phonon process is com- 
parable in magnitude to that from bremsstrahlung from the static lattice. 

The phonon processes, which at high temperatures are primarily Umklapp ones, begin 
to be cut off at a temperature for which the wavenumber of a thermal phonon is comparable 
with the range of wavenumbers for which the electron spectrum is significantly distorted 
from the free-electron form in the vicinity of a band gap. This amounts to the condition 
that the thermal energy must be less than or of order Cg/c times the band gap, where 
is a typical sound speed |T^. Thus phonon- induced bremsstrahlung will be suppressed at 



temperatures which are smaller by a factor of order Cg/c than the temperatures at which 
bremsstrahlung from a static lattice is suppressed. For the lowest reciprocal lattice vectors, 
the corresponding temperature is of order TpZ^^^a, while for the largest reciprocal lattice 
vectors it is of order TpZ~^^^a. At the lowest temperatures phonon-induced bremsstrahlung 
of neutrinos will take place by normal processes, as opposed to Umklapp ones, and the 
emission will vary as T^^ [Q. Thus as the temperature is lowered, phonon-induced processes 
will become more important that bremsstrahlung from the static lattice. However, as one 
can see from Fig. 6, where the rate for emission of neutrino pairs by phonons is shown, 
calculated from Eq.(44), the phonon rate remains below the rate for the static lattice at 
temperatures above about 10® K. We remark that if one were to use an expression more 
accurate than Eq.(44) for the rate of the phonon process, this conclusion would not be 
altered. In a future paper we shall make a more detailed comparison of the two processes 
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V. CONCLUSION 

The most important lesson to be learned from these calculations is that, even though 
the electron-ion interaction has a small effect on the energy per electron, since it is of order 
Z'^^^aHe per electron, it can have a major effect on kinetic process at temperatures low 
compared with band gaps. For the process we have considered in this paper, the effects 
are large because they involve interband transitions. For processes which are mediated by 
intraband transitions the effects will generally be milder because the electron states in the 
vicinity of band gaps will be relatively less important. 

Central to our calculations is the nature of the electronic band structure for ultra- 
relativistic electrons. This is very similar to the non-relativistic case treated in standard 
texts on condensed matter physics, but with the important difference that for the ultra- 
relativistic case the matrix elements of the potential contain the familiar helicity factor 
which cuts down the amplitude for large-angle scattering. This factor has previously been 



considered in the scattering of relativistic electrons by phonons [|Tq] . 

The basic conclusion of our calculations is that the rate of emission of electron pairs by 
electrons is strongly suppressed at temperatures below the band gap. For the inner crust of 
neutron stars, the suppression is substantial at temperatures of order 10^ K, or precisely in 
the temperature range where it was earlier argued that it would be the dominant neutrino 
emission process should the processes in the core be suppressed by superfluidity of the 
neutrons and protons, which are estimated to have transition temperatures that can range 
up to about 10^° K. Of course, estimates of superfluid transition temperatures are uncertain, 
but our results suggest that the range of temperatures where neutrino bremsstrahlung by 
electrons could play a significant role in neutron star cooling are extremely limited. 

There are a number of issues that need to be addressed in order to arrive at better 
estimates of the rate of neutrino bremsstrahlung emission by electrons. A major one is the 
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nature of the collective oscillations of matter in the inner crust which enter the estimates 
of Debye- Waller factors. At the lowest densities, the density of neutrons outside nuclei is 
low and they have little effect on the collective modes, which are just the vibrations of the 
ionic lattice, with a neutralizing background of electrons. At higher densities, the role of 
the neutrons becomes more important, and one must consider a lattice of nuclei coupled 



to a neutron fluid, a problem that has been addressed by Sedrakian |2^. However, the 
problem is not a trivial one, as may be seen from the fact that the most straightforward 
estimate of the effective mass of a nucleus moving through a normal Fermi liquid leads to 



a divergent answer |2l|. The collective mode structure of the liquid-crystal-like phases with 
non-spherical nuclei is a completely open problem, but these modes will come to resemble 
the corresponding giant resonances in nuclear matter as the transition to the uniform phase 
is approached. 
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FIG. 2. The basic bremsstrahlung process. The double Hne is the propagator for a band 
electron. 
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FIG. 5. Energy emission rates, for several values of v\\ = K/{2kF), compared with the corre- 
sponding quantities in the high temperature limit. Solid lines: Form factor Fg = 1 ; Dashed lines: 
Form factor Fg = exp{—q'^rjf/10), vn = 8.7fm. 
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FIG. 6. Solid line: Energy emission rate, Eq.(28), compared with the high temperature limit, 
Eq.(40), as a function of temperature, with form factor Fq = 1. Alternating short and long dashes: 
Same as above, but taking only the contribution of the smallest reciprocal lattice vectors to Eq. (28) . 
Short dashes: Same as solid line, but with form factor Fq = exp(— r^v = 8.7fm. Long 
dashes: The emission rate due to conversion of thermal phonons into neutrino pairs, Eq.(44), as 
compared with Eq.(40). 
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